D-wave overlapping band model for cuprate superconductors 



o 
o 

(N 
oo 



o 
o 
I 

Oh 



O 

o 



> 

m 
m 

o 

O 



X 



Susana Orozco, Rosa Maria Mendez-Moreno, Maria de los Angeles Ortiz, and Gabriela MurguicEl 

Departamento de Fisica, Facultad de Ciencias, 
Universidad Nacional Autonoma de Mexico, 
Apartado Postal 21-092, 04021 Mexico, D. F., Mexico 

Within the BCS framework a multiband model with d-wave symmetry is considered. Generalized 
Fermi surface topologies via band overlapping are introduced. The band overlap scale is of the order 
of the Debye energy. The order parameters and the pairing have d-wave symmetry. Experimental 
values reported for the critical temperatures Tc{x) and the order parameters, Ao(a;), in terms of 
dopping X are used. Numerical results for the coupling and the band overlapping parameters in 
terms of the doping are obtained for the cuprate superconductor La2^xSrxCu04. 
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I. INTRODUCTION 

Measurements of angle-resolved photoemission spec- 
troscopy (ARPES)|lj] and tunneling 3, provide enough 
evidence for the relevant role of phonons in high-Tj, super- 
conductivity (HTSC). Experimental data accumulated so 
far for the high- Tc copper-oxide superconductors have 
given some useful clues to unravel the fundamental in- 
gredients responsible for the high transition temperature 
Tc- However, the underlying physical process remains 
unknown. In this context, it seems crucial to study new 
ideas that use simplified schematic models to isolate the 
mechanism(s) that generate HTSC. 

Pairing symmetry is an important element toward 
understanding the mechanism of high-Tc superconduc- 
tivity. Although early experiments were consistent 
with s-wave pairing symmetry, recent experiments sug- 
gest an anisotropic pairing behavior Q. For many 
cuprate superconductors it is generally accepted that 
the pairing symmetry is d-wave for hole-doped cuprate 
superconductors [1] as for electron doped cupratesQ]. On 
the other hand, recent experiments with Raman scatter- 
ing and ARPESji, 0] have shown that the gap struc- 
ture on high-Tc cuprate superconductors, as a func- 
tion of the angle, is similar to a d-wave gap[l, Q. 
The small but non-vanishing isotope effects in high- 
Tc cuprates have been shown compatible with d-wave 
superconductivity [To| . A phonon-mediated d-wave BCS 
like model has recently beenpresented to describe layered 
cuprated superconductors [1 1| . The last model account 
well for the magnitudes of Tc and the oxygen isotope 
exponent of the superconductor cuprates. Calculations 
with BCS theory and van Hove scenario have also been 
done with d-wave pairing The validity of d-wave 
BCS formalism in high-Tc superconductor cuprates has 
been supported by measurements of transport properties 
and ARPES[l3. 

Numerous indications point to the multiband nature 



of the superconductivity in doped cuprates. The agree- 
ment of the multiband model with experimental findings, 
suggests that a multiband pairing is an essential aspect 
of cuprate superconductivity[l3|. 

First principle calculations show overlapping energy 
bands at the Fermi level [l^. The short coherence length 
observed in high-Tc superconductors, has been related 
to the presence of overlapping energy bands [iTj . A 
simple model with generalized Fermi surface topologies 
via band overlapping has been proposed based on indi- 
rect experimental evidence. That confirms the idea that 
the tendency toward superconductivity can be enhanced 
when the Fermi level lies at or close to the energy of a sin- 
gularity in the density of states (DOS)[l^. This model 
that can be taken as a minimal singularity in the density 
of states and the BCS framework, can lead to higher Tc 
values than those expected from the traditional phonon 
barrier. In our model, the energy band overlapping, mod- 
ifies the DOS near the Fermi level allowing the high Tc 
values observed. A similar effect can be obtained with 
other mechanisms as a van Hove singularity in the den- 
sity of states figf. 

The high-Tc copper-oxide superconductors have a char- 
acteristic layered structure: the Cu02 planes. The 
charge carriers in these materials are confined to the two 
dimensional (2D) Cu02 layers [20l|. These layering struc- 
tures of high-Tc cuprates suggest that two-dimensional 
physics is important for these materials 11 1 . 
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In this work, within the BCS framework, a phonon 
mediated d-wave model is proposed. The gap equation 
(with d-wave symmetry) and two-dimensional general- 
ized Fermi surface topologies via band overlapping are 
used as a model for HTSC. A two overlapping band model 
is considered as a prototype of multiband superconduc- 
tors. For physical consistency, an important requirement 
of the model is that the band overlapping parameter is 
not larger than the cutoff Debye energy, Ed- The model 
here proposed will be used to describe some properties of 
the cuprate superconductor La2~xSrxCuOA in terms of 
the doping and the parameters of the model. 
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II. THE MODEL 



We begin with the famous gap equation 



in the weak couphng hmit, with V{k,k') the pairing 
interaction, ks is the Boltzman constant, and — 
el + Af., where e/c = ?i^fc^/2m are the self-consistent 
single-particle energies. 

For the electron-phonon interaction, we have consid- 
ered, with Vb a constant, V{k, k') = Voip{k)'tp{k') when 
|efe| and \ek'\ < Ed — ksTo and elsewhere. As usual 
the attractive BCS interaction is nonzero only for unoc- 
cupied orbitals in the neighborhood of the Fermi level 
Ep. In the last equation, tp{k) = cos(2 0^) for d^2_y2 
pairing. Here 0^ = t&'D~^ {ky / kx) is the angular direction 
of the momentum in the ab plane. The superconducting 
order parameter, A(fc) = A(r) ijj{k) if |efe| < En and 
elsewhere. 

With these considerations we propose a generalized 
Fermi surface. The generalized Fermi sea proposed con- 
sists of two overlapping bands. As a particular distri- 
bution with anomalous occupancy in momentum space 
the following form for the generalized Fermi sea has been 
considered 

rife = e(7fcF - A:) + Qi^kp - k)Q{k - pkp), (2) 

with kp the Fermi momentum and < /? < 7 < 1. 
In order to keep the average number of electron states 
constant, the parameters are related in the 2D system by 
the equation 



272 - = 1, 



(3) 



then only one of the relevant parameters is independent. 
The distribution in momentum induces one in energy, 
Efj < Ey where E/s = (3'^Ep and E^ — j^Ep . We require 
that the band overlapping be of the order or smaller than 
the cutoff (Debye) energy, which means (1-7^)^;;. < Ed. 
The last expression can be written as 



(1 - j^)Ep = t^Ed, 



(4) 



where 77 is in the range < 77 < Ep /{2Ed)- Equations 
([5]) and (HI together will give the minimum 7^ value con- 
sistent with our model. 

In the last framework the summation in Eq. ^ is 
changed to an integration which is done over the (sym- 
metric) generalized Fermi surface defined above. One 
gets 
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In this equation = (efc - Ep^ + A(T)2 cos2(2 0), 
the coupling parameter is A = VqD{E), with D{E) the 
electronic density of states, which will be taken as a 
constant for the 2D system in the integration range. 
Ep = ^^n2Dj with n2D the carriers density per CUO2 
layer. The two integrals correspond to the bands pro- 
posed by Eq. (g]). 

The integration over the surface at E-y in the first 
band, is restricted to states in the interval Ej — Ed < 
Ek < Ej + Ed- In the second band, in order to con- 
serve the particle number, the integration is restricted 
to the interval Ep < Ek < Ep, if E^ + Ed > Ep, 
with Ep = (2 7^ - l)Ep, according to Eq. ^ 
in our model. While Ep ~ E^ < Ed, implies that 
the energy difference between the anomalously occupied 
states must be provided by the material itself. Finally 
A(T) ipik) = A(r) cos(2 0) at the two bands. 

The critical temperature is introduced via the Eq. ([S]) 
at T = Tc, where the gap becomes A.{Tc) = 0. At this 
temperature Eq. ([5]) is reduced to 
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tanh 
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Efc — Ef 



dek 



2kBTc ) ek - Ep 
ek- Ep\ dek 



(6) 



2kBTc I ek - Ep 



which will be numerically evaluated. The last equation 
relates Tc to the coupling constant A and to the anoma- 
lous occupancy parameter 7^. This relationship deter- 
mines the 7^ values which reproduces the critical tem- 
perature of several cuprates in the weak coupling region. 

At T = OK, Eq. (O will also be evaluated and 7^ values 
consistent with the numerical results of Eq. ([6]) will be 
obtained: 
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(7) 



sinh" 



Ao |cos(2 0)1 



sinh- (1 - 



(5) 



Ao I cos (2 </))! 



where A(0) = Aq. 

The model presented in this section can be used to de- 
scribe high-Tc cuprate superconductors, the band over- 
lapping 1 — 7^ and relevant parameters are determined. 
In any case a specific material must be selected to in- 
troduce the available experimental data. Ranges for the 
coupling parameter A in the weak coupling region, and 
the overlapping parameter 7^, consistent with the model 
and the experimental data, can be obtained for each ma- 
terial. The relationship between the characteristic pa- 
rameters will be obtained for La-based compounds at 
several doping concentrations x, ranging from the under- 
doped to the overdoped regime. Different values of the 
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coupling constant and the overlapping parameter consis- 
tent with the model, are obtained using the experimental 
values of Aq and Tc. 

The single layer cuprate superconductor 
La2-xSrxCuOi [La — 214) has one of the simplest 
crystal structures among the high-Tc superconductors. 
This fact makes this cuprate very attractive for both 
theoretical and experimental studies. High quality single 
crystals of this material are available with several doping 
concentrations which are required for experimental 
studies. Even the determination of charge carrier 
concentration in the cuprate superconductors is quite 
difficult, the La — 214 is a system where the carrier 
concentration is nearly unambiguously determined. For 
this material, the hole concentration for Cu02 plane, 
n2D, is equal to the x value, i.e. to the Sr concen- 
tration, as long as the oxygen is stoichiometric [2li [23 |. 
Additionally, there are reliable data for the and 
the superconducting gap Aq for several samples in the 
superconducting region. 
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FIG. 1: The coupling parameter A in terms of the overlapping 
parameter 7^^, with optimal doping x — 0.16. The T = OK 
(broken curve) and T = Tc = 40K (continuous curve) for 
La2-xSrxCu04,, are shown. The horizontal line at A = 0.68 
shows the maximum A selected for T = Tc with 7^ = 0.95. 



III. RESULTS AND DISCUSSION 

In order to get numerical results, with our overlap- 
ping band model with d-wave symmetry, the cuprate 
La2-xSrxCuOi was selected. The values for Aq are 
taken in the interval 2 < Aq < 12 meV which includes 
experimental results [22|. The behavior of A as function 
of X and 7^ at T = Tc is obtained from Eq. ([6]) ; and A as 
function of Ao, x and 7^ at T = OK is given by Eq. ([7]). 
To have coupled solutions of these equations the same A 
value for T — Tc and T — OK is proposed. These so- 
lutions correspond to different overlap values 1 — 7^, at 
each equation. With this model and s-wave symmetry, 
the band overlapping 1 — 7^ was higher at T = OK than 
at T = Tc[2^. We consider the same behavior with d- 
wave symmetry. The maximum Tc for cuprate supercon- 
ductors is obtained at optimal doping. With the model 
A(a;) values are obtained, including at optimal doping 

X{Xop)M- 

In Fig. [TJ values of the coupling parameter A in terms 
of the overlapping parameter 7^ are shown in the weak 
coupling re gion . The experimental results of Tc and Aq 
from Refs. [2^ and [l^l were introduced. The curves at 
T = OK (broken curve) and at T = Tc = 40K (con- 
tinuous curve) for La2-xSrxCu04, with optimal doping 
Xop = 0.16 are shown. The minimum 7^ value of 0.55 
was taken to be consistent with the model. In the whole 
range reported for the band overlapping, the coupling 
parameter required at each 7^ is larger for T = OK than 
for T = Tc. In order to use the same A for T — Tc and 
T = OK, the A values must be restricted i.e., the A value 
at each 7^ must be larger than A,„i„ — 0.57 at T — OK. 



In the region 7^ > 0.7 with a constant A value, a larger 
band overlapping 1 — 7^ is obtained for T = OK than for 



T = Tc in agreement with our assumption. For instance, 
the maximum A for T = Tc with 7^ = 0.95, is shown by 
the horizontal line at A = 0.68, and the intersection of 
this line and the T — OK curve is at 7^ — 0.94. The same 
restrictions over A are considered at any other doping in 
the superconducting phase. However, for any x 7^ Xop, 
the A value must be smaller than A = 0.68. 

In Fig. [2] the results for optimal doping Xop, are com- 
pared with the underdoped x — 0.13 and the overdoped 
X = 0.2 cases. The experimental values of Aq and Tc 
for each doping, were introduced. The continuous curves 
correspond to Xop, the small dashed curves show the un- 
derdoped behavior and the large dashed ones the over- 
doped results. In the optimal doped and underdoped 
cases, the T = OK curves are above the corresponding 
T — Tc ones. In the overdoped case, the behavior is 
different i.e., the Tc curve is above the T = OK one. 



In the three cases, the values of the coupling param- 
eter are in the weak coupling region for the 7^ values 
which satisfy the conditions of our model. All the 7^ 
values which satisfy the A restrictions are allowed. How- 
ever, as an example, we have selected extreme A values 
in the three cases. The three horizontal lines show these 
A values. 

As in Fig. [1] the maximum A value selected at optimal 
doping is A = 0.68. In the underdoped case A — 0.65 is 
selected. This value corresponds to the overlapping pa- 
rameter 7^ = 0.621, the minimum of the T = OK curve, 
and 7^ = 0.941 at the T — Tc curve. In the overdoped 
case, the selected A value is 0.51, i.e. the minimum of the 
curve T = Tc. With this A, the overlapping parameters 
are 7^ = 0.599 for T = OK and 7^ = 0.76 for T = Tc. 

With numerical solutions of Eq. ([7]) we may obtain the 
gap Aq in terms of the parameters of our model. The 
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FIG. 2: The results for optimal doping (continuous curves) 
are compared with the underdoped (small dashed curves) x — 
0.13 and the overdoped x = 0.2 (large dashed curves) cases. 
The three horizontal lines show the selected extreme A values: 
A — 0.68 at optimal doping, A — 0.65 in the underdoped case, 
and A = 0.51 in the overdoped case. 

underdoped material is considered in Fig. [3] because the 
advantage of our model is easily shown. The gap Aq is 
shown in terms of the coupling parameter A. The gap 
Aq always increases with the coupling parameter A. The 
curves are drawn for 7^ — 0.621, 0.5 and 0.8 from up to 
down respectively. For this sample, with 7^ = 0.621, we 
obtain the minimum A value for any Aq and for any A 
the maximum Ao value. 
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FIG. 3: The gap Ao obtained in terms of the coupling 
parameter A for the underdoped sample. The curves are 
drawn for 7^ — 0.621,0.5 and 0.8 from up to down respec- 
tively. The continuous horizontal line shows the experimental 
Ao = 10.85meV value. The large dashed horizontal line shows 
the d-wave mean-field approximation Aa/f = 6.68meV result. 



The continuous horizontal line shows the experimen- 
tal Aq = 10.85meV value. The large dashed hori- 
zontal line shows the d-wave mean-field approximation 



Amf = 6.68meV result [25"], where the same d-wave sym- 
metry was considered. However, introduction of the band 
overlapping allows to reproduce the experimental result 
with all the 7^ values in the range considered. The band 
overlapping model also allows higher Aq values for the 
underdoped system and lower Aq for the overdoped one, 
than the I^mf = 2.14:5kBTc. 

In Fig. m the behavior between Aq and A for optimal 
doping is compared with the underdoped and the over- 
doped cases. The 7^ values introduced are those selected 
in Fig.[2]for T — OK. The horizontal lines are the A values 
also selected in Fig. [21 All the continuous curves corre- 
spond to optimal doping. The large and small dashed 
curves correspond to the overdoped and the underdoped 
systems respectively. The curves show the interesting 
relationship between these parameters. As for optimal 
doping, the coupling parameter increases with Aq for any 
doping. The vertical lines are the experimental Aq val- 
ues. It is possible to reproduce the experimental Aq in 
the range 0.13 < x < 0.2. The band overlapping in- 
troduced in this model allows the reproduction of the 
behavior of Aq with doping. 
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FIG. 4: Ao as function of A for optimal doping, underdoped, 
and the overdoped cases. The horizontal lines are the A values 
selected in Fig. [2] The vertical lines show the experimental 
Ao values. The continuous lines correspond to a; = 0.16, the 
small dashed ones to a; = 0.13, and the large dashed ones to 
X = 0.2. 



In conclusion, we presented an overlapping band model 
with d-wave symmetry, to describe high- Tc cuprate su- 
perconductors, within the BCS framework. We have used 
a model with anomalous Fermi Occupancy and d-wave 
pairing in the 2D fermion gas. The anomaly is introduced 
via a generalized Fermi surface with two bands as a pro- 
totype of bands overlapping. We report the behavior of 
the coupling parameter A as function of the gap Aq and 
the overlapping parameter 7^, for different doping sam- 
ples. The A values consistent with the model are in the 
weak coupling region. The behavior of Aq as function 
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of A shows that for several band overlapping parameters 
it is possible to reproduce the experimental Aq values 
near the optimal doping, for the cuprate La2-xSrxCu04. 
The band overlapping allows the improvement of the re- 
sults obtained with a d-wave mean-field approximation, 
in a scheme in which the electron-phonon interaction is 
the relevant high-Tg mechanism. The energy scale of the 



anomaly (1 — 7^)ii^i^ is of the order of the Debye energy. 
The Debye energy is then the overall scale that deter- 
mines the highest and gives credibility to the model 
because it requires an energy scale accessible to the lat- 
tice. The enhancing of the DOS with this model simu- 
lates quite well intermediate and strong coupling correc- 
tions to the BCS framework. 
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